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Abstract 

Canine distemper virus (CDl/9-induced demyelinating leukoencephalitis in dogs {Canis familiaris) is suggested to represent a 
naturally occurring translational model for subacute sclerosing panencephalitis and multiple sclerosis in humans. The aim of 
this study was a hypothesis-free microarray analysis of the transcriptional changes within cerebellar specimens of five cases 
of acute, six cases of subacute demyelinating, and three cases of chronic demyelinating and inflammatory CDV 
leukoencephalitis as compared to twelve non-infected control dogs. Frozen cerebellar specimens were used for analysis of 
histopathological changes including demyelination, transcriptional changes employing microarrays, and presence of CDV 
nucleoprotein RNA and protein using microarrays, RT-qPCR and immunohistochemistry. Microarray analysis revealed 780 
differentially expressed probe sets. The dominating change was an up-regulation of genes related to the innate and the 
humoral immune response, and less distinct the cytotoxic T-cell-mediated immune response in all subtypes of CDV 
leukoencephalitis as compared to controls. Multiple myelin genes including myelin basic protein and proteolipid protein 
displayed a selective down-regulation in subacute CDV leukoencephalitis, suggestive of an oligodendrocyte dystrophy. In 
contrast, a marked up-regulation of multiple immunoglobulin-like expressed sequence tags and the delta polypeptide of the 
CD3 antigen was observed in chronic CDV leukoencephalitis, in agreement with the hypothesis of an immune-mediated 
demyelination in the late inflammatory phase of the disease. Analysis of pathways intimately linked to demyelination as 
determined by morphometry employing correlation-based Gene Set Enrichment Analysis highlighted the pathomechanistic 
importance of up-regulated genes comprised by the gene ontology terms "viral replication" and "humoral immune 
response" as well as down-regulated genes functionally related to "metabolite and energy generation". 
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Introduction 

Canine distemper virus (CDV) is a morbillivirus of the family 
Paramyxoviridae and the etiologic agent of distemper in dogs {Canis 
familiaris). Distemper is a naturally and worldwide-occuring, highly 
contagious systemic virus infection of dogs as well as various other 
carnivores [1]. The mortality rate can be up to 50% but is often 
much lower [2]. The route of infection is mainly oronasally 
followed by a first viremic spread of CDV to lymphoid tissues 
resulting in lymphopenia and immunosuppression [3]. Subse- 
quently, depending on the virulence of the respective CDV strain, 
the host's age and immune status, dogs may either initiate a robust 
humoral immune response and recover or may develop a second 
viremic spread to multiple organs and tissues [4,5,6]. These dogs 
show various pathologic alterations including thymic atrophy, 
conjunctivitis, rhinitis, interstitial pneumonia, encephalitis, gastro- 
enteritis, pustular dermatitis and hyperkeratosis of the foot pads 



and nose [3,7,8,9]. Dogs may either succumb to CDV infection 
during this phase or recover. However, some dogs may develop a 
chronically persistent or late-onset encephalitis [10,11]. 

CDV-induced encephalitis frequently manifests as an acute, 
subacute or chronic leukoencephalitis, and rarely as a polioen- 
cephalitis [12,13]. CDV leukoencephalitis lesions are found mainly 
periventricularly in the cerebellum, optic tract, medulla oblongata, 
cerebral peduncle, spinal cord, and less frequently in the 
subcortical white matter, fornix, capsula interna and corpus 
callosum [12,14]. Based on histopathological findings in experi- 
mentally-induced CDV leukoencephalitis in dogs, central nervous 
system (CNS) lesions can be classified into distinct subtypes 
presumably developing in sequential order [14,15]. Briefly, the 
hallmark of acute lesions are focal accumulations of CDV antigen- 
positive glial cells, eventually accompanied by vacuolation of the 
white matter, and mild gliosis [14,16,17]. This type of lesions has 
been observed 16-24 days after experimental CDV-infection 
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Table 1. Experimental design and signalement of dogs with CDV leukoencephalitis and controls. 



Group ID 


Animal ID 


Breed 


Gender 


Age (month) 


Body weight (kg) 




V723-06 


Beagle 


m 


6 


8.25 




V725-06 


Beagle 


f 


6 


10.80 




V776-06 


Beagle 


f 


5 


6.50 




V777-06 


Beagle 


m 


5 


7.45 




V239-05 


Beagle 


m 


19 


n.d. 




V255-05 


Beagle 


m 


20 


n.d. 




V256-05 


Beagle 


f 


17 


n.d. 




V259-05 


Beagle 


f 


17 


n.d. 




S 1486-98 


Cavalier King Charles Spaniel 


f 


8 


3.00 




S1 497-98 


Rhodesian Ridgeback 


m 


8 


32.00 




S521-99 


Giant Schnauzer 


f 


5 


16.00 




SI 863-98 


Yorkshire Terrier 


f 


3 


0.75 


2 


S 104-95 


Mixed breed 


m 


5 


5.50 


2 


S2548-90 


Bobtail 


m 


6 


n.d. 


2 


S842-04 


Mixed breed 


f 


4 


5.00 


2 


S997-03 


Jack-Russell-Terrier 


m 


5 


3.00 


2 


S580-05 


Mixed breed 


m 


3 


2.10 


3 


S1 358-90 


Terrier 


f 


4 


4.00 


3 


S98-91 


Mixed breed 


f 


6 


n.d. 


3 


S2015-91 


German sheperd 


m 


4 


n.d. 


3 


S1 03-95 


Dachshund 


f 


12 


5.00 


3 


S58-05 


Husky 


f 


7 


12.80 


3 


S954-05 


Mixed breed 


m 


6 


13.00 


4 


S2936-04 


Jack-Russell-Terrier 


f 


3 


3.50 


4 


S1510-90 


Mixed breed 


m 


47 


n.d. 


4 


S2290-90 


German sheperd 


m 


10 


n.d. 



Group ID: 1 = control dogs, 2 = acute CDV leukoencephalitis, 3 = subacute CDV leukoencephalitis, 4 = chronic CDV leukoencephalitis; f = female; ID = identifier; m = male; 

n.d. = not determined. 

doi:1 0.1 371 /journal.pone.009591 7.t001 



[15,16,18,19]. Subacute lesions with demyelination but without 
inflammation may occur 24—32 days post infection 
[15,16,18,19,20]. Subacute to chronic lesions with demyelination 
and inflammation are characterized by reduced numbers of CDV 
antigen-positive cells and distinct perivascular lymphohistioplas- 
macytic cuffs of several layers of thickness [14,16,17]. Following 
experimental infections, this type of lesions occurred after a 
minimum of 29-63 days post infection [5,15,16,18,19,20], Overall 
the correlation between morphological changes and assumed time 
point of infection remains a vague supposition in cases of 
spontaneously occurring CDV leukoencephalitis, since multiple 
factors including virus strain, immune status and age of the host as 
well as various environmental parameters influence virus spread 
and disease progression [4,5,6]. 

CDV leukoencephalitis in dogs has been suggested as a 
spontaneously-occuring viral model to study the pathogenesis of 
demyelination in subacute sclerosing panencephalitis and multiple 
sclerosis in humans [21,22]. Suggested mechanisms mediating 
primary demyelination in CDV leukoencephalitis include virus- 
induced damage of oligodendrocytes [19,23], bystander damage 
due to activated microglia/macrophages releasing myelinotoxic 
reactive oxygen species and proteolytic enzymes such as matrix 
metalloproteinases [24,25,26], activated astrocytes and microglia/ 



macrophages releasing tumor necrosis factor-oc (TNFot) [27], and 
humoral and cell-mediated autoimmunity [28,29,30,31]. 

In order to shed more light on the molecular pathogenesis of 
CDV leukoencephalitis, the aims of the present investigation were 
1.) to perform a classic assumption-free microarray analysis of 
transcriptional changes in CDV leukoencephalitis, and 2.) to 
identify pathways intimately associated with demyelination by 
combining findings obtained by histological analysis of brain 
sections of different subtypes of CDV leukoencephalitis with 
respective microarray data employing correlation-based Gene Set 
Enrichment Analysis. 

Materials and Methods 

Ethics Statement 

This study was conducted in accordance with the German 
Animal Welfare Act. The authors confirm that no animals were 
infected or sacrificed for the purpose of this retrospective 
pathological case-control study. This study is not an animal 
experiment since all animals were dead at the time of submission 
for necropsy in order to investigate the causes of death and disease. 
All tissues used in this study were collected by one of the authors 
(WB) during his work at the diagnostic pathology services of the 
Department of Pathology, University of Veterinary Medicine 
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Figure 1. Pathohistological changes characteristic for the different subtypes of CDV leukoencephalitis. (A) The cerebella of the non- 
infected control dogs (group 1) displayed no histological alterations. (B) The cerebella of dogs affected by acute CDV leukoencephalitis (group 2) 
showed focal astro- and microgliosis (arrow) and occasionally few vacuolated myelin sheaths. (C) The cerebella of dogs affected by subacute CDV 
leukoencephalitis with demyelination but without inflammation (group 3) exhibited focally demyelinated white matter (asterix) combined with astro- 
and microgliosis (arrow). (D) The cerebella of dogs affected by chronic CDV leukoencephalitis with demyelination and with inflammation (group 4) 
displayed focally demyelinated white matter (asterix), combined with astro- and microgliosis (arrow) as well as perivascular inflammatory infiltrates 
(arrowhead). Luxol fast blue-cresyl violet. Scale bars = 200 urn. 
doi:10.1371/journal.pone.0095917.g001 



Hannover, and the Institute of Veterinary Pathology, Justus- 
Liebig-University Giessen, and most animals were used in previous 
publications [32,33]. All dog owners provided written consent for 
the dogs' tissues to be collected and used for research purposes. 

Experimental Design 

Age, breed and gender of the dogs are presented in Table 1 . 
Twelve control animals (group 1) covering a comparable range of 
variation with respect to gender, age and breed as CDV-infected 
dogs, and lacking morphological evidence of CNS disease were 
selected from the departmental archives. The 14 CDV-infected 
dogs suffered from spontaneously-occurring and immunohistolog- 
ically-confirmed CDV leukoencephalitis. For the current study, 
only specimens of cerebellum and adjacent brain stem unequiv- 
ocally displaying focal or multifocal lesions of a single subtype of 
CDV leukoencephalitis within each individual dog were selected. 
Based on histopathologic al findings, dogs were classified into acute 
CDV leukoencephalitis (group 2, N = 5), subacute CDV leukoen- 
cephalitis with demyelination but without inflammation (group 3, 
N = 6), and (subacute to) chronic CDV leukoencephalitis with 



demyelination and inflammation (group 4, N = 3) according to 
previously described criteria [17,33]. 

Histology and Immunohistology 

Cerebellar specimens were snap-frozen in Tissue-Tek O.C.T. 
compound (Sakura Finetek Europe B.V. Zoeterwoude, Nether- 
lands) and ~4 |0.m thick sections, contiguous to the tissue used for 
RNA isolation, were cut on a cryotome, mounted on SuperFrost 
Plus microscope slides (Menzel-Glaser, Braunschweig, Germany), 
and stained with hematoxylin and eosin (HE), and Luxol fast blue- 
cresyl violet. Immunohistological examination of serial sections 
was performed with the avidin-biotin-peroxidase complex (ABC) 
method (Vector PK-6100) and 3,3'-diaminobenzidine-tetrahy- 
drochloride as chromogen. A polyclonal rabbit antibody specific 
for the CDV nucleoprotein (diluted 1:2000, #25, kindly provided 
by C. Orvell, The National Bacteriological Laboratory, Stock- 
holm, Sweden) was used to confirm CDV infection [33]. 
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780 differentially expressed probe sets 




Figure 2. Heatmap displaying the expression profile of the 780 DEPs. Hierarchical cluster analysis of the 780 DEPs was performed employing 
Euclidean distance and complete linkage to reveal similar expression patterns. Each column represents one of the 780 DEPs and each row one of the 
26 biological replicates (cerebellar specimens of individual dogs) sorted according to the histologically defined subtypes of CDV leukoencephalitis: 
group 1 = controls; group 2 = acute CDV leukoencephalitis; group 3 = subacute CDV leukoencephalitis; group 4 = chronic CDV leukoencephalitis. The 
heatmap displays the log 2 -transformed individual fold changes relative to the mean expression of the controls indicated by a color scale ranging from 
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-2 (4-fold down-regulation) in green to 2 (4-fold up-regulation) in red. Accordingly, the DEPs are subdivided into four clusters with distinct expression 
profiles (yellow, brown, orange and green bars below the heatmap). Gene ontology analysis of the orthologous mouse genes followed by functional 
annotation clustering revealed distinct biological modules within the DEPs of cluster 1, 3, and 4. * = important associated genes and expressed 
sequence tags are reported for the DEPs of cluster 2 since automated gene ontology analysis revealed no significant results due to poor annotation 
quality. ES = enrichment score. 
doi:1 0.1 371 /journal.pone.009591 7.g002 



Table 2. Transcriptional differences between controls (group 1) and subtypes of CDV leukoencephalitis (group 2-4). 





Pairwlse 

post-hoc 
tests 1 


DEPs 
total 


DEPs 
up/down 


DEGs* up/ 
down 


Biological modules' 


Group 2 
versus 1 


423 


up 398 


231 


positive regulation of immune response (ES 19.71); 
immunoglobulin mediated immune response (ES 16.82); 
complement activation, classical pathway (ES 9.88); 
positive regulation of lymphocyte activation (ES 5.04); 
positive regulation of T-cell differentiation (ES 3.83); 
positive regulation of T-cell mediated cytotoxicity (ES 2.65); 
positive regulation of B-cell proliferation (ES 2.51) 






down 25 


16 


regulation of apoptosis (ES 2.10); positive regulation 
of apoptosis (ES 1 82)' central nervous system 
development (ES 1.30); positive regulation of transcription 
from RNA polymerase II promoter (ES 1.23) 


Vj 1 U U U ' 

versus 1 


688 


up 549 


31 1 


r»nciti\/fi rf^ni i ^tinn nf immi i n rocnnncp 1 ft QAV 

ULolllVtr I *r UJ U I cH J KJ I I KJ\ 1 1 I 1 1 I I U I I ^ [ L. J ^JKJ 1 1 {I- -J 1 O ■ z/\ ) , 

immunoglobulin mediated immune response (ES 14.98); 
complement activation, classical pathway (ES 8.63); 
positive regulation of lymphocyte activation (ES 5.21); 
positive regulation of T-cell differentiation (3.95); 
positive regulation of T-cell mediated cytotoxicity (ES 2.42) 






down 139 


75 


intermediate filament bundle assembly (ES 2.82); 
fatty acid biosynthetic process (ES 1.28) 


Group 4 
versus 1 


525 


up 486 


280 


immunoglobulin mediated immune response (ES 12.21); 
complement activation, classical pathway (ES 7.03); 
positive regulation of lymphocyte activation (ES 5.58); 
positive regulation of T-cell differentiation (ES 4.42); 
positive thymic T-cell selection (ES 3.15); 
positive regulation of T-cell mediated cytotoxicity 
(ES 2.65); positive regulation of B-cell proliferation 
(ES 2.45); positive regulation of NF-kappaB 
transcription factor activity (ES 2.25) 






down 39 


28 


central nervous system development (ES 1.81}; 
positive regulation of transcription from RNA 
polymerase II promoter (ES 0.82} 


uroup 5 
versus 2 


1 06 




44 


protein amino acid N-linked glycosylation (ES 1.81); 
negative regulation of apoptosis (ES 1.57); 
regulation of lymphocyte proliferation (ES 1.49) 






down 36 


21 


intermediate filament bundle assembly (ES 3.33); 
nervous system development (ES 1.83); 
regulation of neuron differentiation (ES 1.70) 


Group 4 
versus 2 


89 


up 77 


47 


anatomical structure development (ES 0.47); 
response to biotic stimulus (ES 0.43); regulation 
of immune response (ES 0.43); protein metabolic 
process (ES 0.17) 






down 12 


8 


immune response (ES 1.31); regulation of 
biological process (ES 0.26); cellular 
metabolic process (ES 0.04) 


Group 4 
versus 3 


229 


up 100 


54 


intermediate filament bundle assembly (ES 2.94); 
regulation of neuron differentiation (ES 1.27) 






down 129 


75 


complement activation, classical pathway (ES 2.32); 
protein amino acid N-linked glycosylation (ES 1.51) 


1 Multigrou 


p test followed by 


independent pairwise 


post-hoc comparisons employing LIMMA with a false discovery rate (Benjamini and Hochberg) filter of q<0.05, and a 



fold change filter of >2.0 or < -2.0. 

"'"Gene as defined by the DAVID knowledgebase. Multiple probe sets may match to a single gene, and probe sets may be unassigned (expressed sequence tags). 

§ Biological modules associated with the DEG as detected by the DAVID functional annotation clustering algorithm using orthologous mouse genes. 

DEGs = differentially expressed genes; DEPs - differentially expressed probe sets; ES = enrichment score; group 1 = control dogs, group 2 = acute CDV leukoencephalitis, 

group 3 -subacute CDV leukoencephalitis, group 4 = chronic CDV leukoencephalitis. 

doi:1 0.1 371 /journal.pone.009591 7.t002 
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defense reponse (ES 8.56) 
response to other organism (ES3.13) 
immunoglobulin mediated immune response (ES1.75) 
induction of apoptosis (ES 1.34) 



cellular protein localization (ES 0.72) 



group 2 versus 1 



positive regulation of immune 
response (ES 18.22) 

immunoglobulin mediated 
immune resonse (ES 13.84) 

complement activation, classical 
pathway (ES 7.83) 

positive regulation of lymphocte 
activation (ES 5.71) 

positive regulation of T-cell 
differentiation (ES 3.96) 

positive regulation of B-cell 
proliferation (ES 2.77) 

lymphocyte differentiation (ES 2.64) 

positive regulation of T-cell 
mediated cytotoxicity (ES 2.63) 




negative regulation of cellular 
biosynthetic process (ES 1.38) 

apoptosis (ES 1.32) 

nduction of apoptosis (ES 1.31) 

regulation of apoptosis (ES 1.18) 

immune system 
development (ES 1.11) 



group 4 versus 1 



regulation of immune response (ES 1.39) 



B 



up-regulated genes 



group 2 versus 1 group 3 versus 1 



central nervous system 
development (ES 1.49) 

positive regulation of 
transcription from RNA 
polymerase II promoter (ES 1.48) 

regulation of transcription, 
DNA dependent (ES 1.36) 

cell differentiation (ES 1.20) 




intermediate filament bundle 
assembly (ES 2.94) 

regulation of neuron 
differentiation (ES 1.50) 



transport (ES 1.17) 



group 4 versus 1 



down-regulated genes 

Figure 3. Venn diagrams comparing the DEGs within the a priori defined subtypes of CDV leukoencephalitis. The three subtypes of 
CDV leukoencephalitis (group 2-4) are compared for shared and unique DEGs as well as associated biological modules as detected by gene ontology 
analysis and functional annotation clustering of the orthologous mouse genes. (A) The Venn diagram displays the numbers and intersections of the 
up-regulated DEGs in the three subtypes of CDV leukoencephalitis (group 2-4) as compared to controls (group 1 ), respectively. (B) The Venn diagram 
displays the numbers and intersections of the down-regulated DEGs in the three subtypes of CDV leukoencephalitis (group 2-4) as compared to 
controls (group 1), respectively. ES = enrichment score; group 1= controls; group 2 = acute CDV leukoencephalitis; group 3 = subacute CDV 
leukoencephalitis; group 4 = chronic CDV leukoencephalitis. 
doi:10.1371/journal.pone.0095917.g003 



RNA Isolation, Microarray Hybridization and Low-level 
Analysis 

Total RNA was isolated from the frozen cerebellar specimens 
using the RNeasy Tipid Tissue Mini Kit (Qiagen, Hilden, 
Germany). The quality and integrity of the RNA was controlled 



using the Agilent 6000 RNA Nano kit and an Agilent Bioanalyzer 
2100 (Agilent, Boblingen, Germany). 200 ng of each RNA was 
amplified and biotin-labeled employing the 3TVT express kit 
(Affymetrix, Santa Clara, USA), and hybridized to GeneChip 
canine genome 2.0 arrays (Affymetrix, Santa Clara, USA) in a 
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Figure 4. Multifocal demyelination is the hallmark of subacute 
and chronic CDV leukoencephalitis. The dot diagram shows an 
increased percentage of Luxol fast blue-negative white matter area 
(myelin loss) in subacute CDV leukoencephalitis (group 3) and chronic 
CDV leukoencephalitis (group 4) as compared to controls (group 1) and 
acute CDV leukoencephalitis (group 2). Each point represents the 
percentage of Luxol fast blue-negative area relative to total white 
matter area within the cerebellar specimen of each individual dog. The 
median of each group is represented by a horizontal line. Significant 
differences between the groups as revealed by the Kruskal-Wallis test 
with independent pairwise post-hoc Mann-Whitney U-tests are marked 
as follows: *p<0.05. 
doi:1 0.1 371 /journal.pone.009591 7.g004 

rotating hybridization oven at 45°C for 16 hrs. Afterwards the 
arrays were washed and stained with a solution containing R- 
phycoerythrin-streptavidin employing the Affymetrix GeneChip 
Fluidics Station 450 (Affymetrix, Santa Clara, USA). Scanning was 
performed with an Affymetrix GeneChip Scanner 3000 (Affyme- 
trix, Santa Clara, USA). Background adjustment, quantile 
normalization and probe set summarization was performed using 
the GC-RMA algorithm (Bioconductor genua for R package, 
Version 2.3) [34]. Principal components analysis was performed as 
a quality check for oudying samples employing Babelomics [35]. 
MIAME compliant data sets are deposited in the ArrayExpress 



database (accession number: E-MEXP-3917; http://www.ebi.ac. 
uk/arrayexpress). The mean expression values of the spiked-in 
hybridization control RNAs (1.5 pM Escherichia coli biotin synthase 
[bioB], 5 pM Escherichia coli biotin biosynthesis protein [bioC], 25 pM 
Escherichia coli dethiobiotin synthetase [bioD], 100 pM bacteriophage PI 
cyclization recombinase [ere], Affymetrix, Santa Clara, USA) were used 
as comparative markers for no to low, moderate, high, or highest 
expression, respectively. 

Differentially Expressed Probe Sets 

Differentially expressed probe sets (DEPs) were detected 
employing a one-factorial test using the Linear Models for 
Microarray Data (LIMMA) algorithm with a maximal false 
discovery rate (FDR) of 5% (q£0.05) according to Benjamini 
and Hochberg [35,36], followed by post-hoc pairwise comparisons 
of group 1-4 combining a statistical significance filter (LIMMA, 
q<0.05) and a fold change filter (fold change >2.0 or < -2.0). 
The fold change was calculated as the ratio of the inverse- 
transformed arithmetic means of the logj-transformed expression 
values. Down-regulations are shown as negative reciprocal values. 
For the hierarchical clusters we used TM4 Multi Experiment 
Viewer (MeV) with the individual fold change of each animal 
relative to the mean of all control dogs, Euclidean distance, and 
complete linkage [37]. Gene lists were compared and analyzed for 
intersections employing Venn diagrams (Oliveros, J.C. VENNY. 
An interactive tool for comparing lists with Venn Diagrams. 
http://bioinfogp.cnb.csic.es/tools/venny/index.html.). 

Annotation and Gene Ontology Information 

Probe sets were annotated with canine gene symbols and gene 
names according to the Affymetrix annotation file (release 33; 29. 
October 2012). Due to non-perfect accordance of the multiple co- 
existing genomic identifier (ID) and nomenclature systems (official 
gene symbol, Unigene ID, Entrez Gene ID, etc.) genes were 
defined according to the integrative DAVID knowledgebase 
(DAVID ID) throughout this text [38]. Orthologous mouse gene 
symbols were retrieved employing MADGene [39]. Lists of DEPs 
were consolidated to lists of differentially expressed genes (DEGs) 



Table 4. Gene ontology terms correlated to demyelination as revealed by Gene Set Enrichment Analysis. 





Gene ontology terms 
(biological process) 


Size 
(N) 


Gene symbols 


Normalized 
enrichment 
score 


P- 

value 


Gene ontology terms positively correlated with Luxol fast b\ue -negative white matter area (myelin-loss) 


Viral reproduction 


5 


CALCOC02, CD81, TOP2A, CCL2, HBXIP 


1.71 


0.0099 


Viral infectious cycle 


4 


CD81, TOP2A, CCL2, HBXIP 


1.65 


0.0224 


Viral reproductive process 


4 


CD81, TOP2A, CCL2, HBXIP 


1.65 


0.0224 


Mitotic cell cycle 


9 


TPX2, KIF23, CDKN1A, BTG3, SMC4, NEK6, CDKN2C, AFAP1L2, CDK2API 


1.62 


0.0358 


Humoral immune response 


4 


BLNK, PSMB10, CCL2, LY86 


1.58 


0.0106 


Regulation of protein metabolic process 


7 


TIM PI, IRAK3, CD8I, TLP.3, DNAJCI, CCND2, TLR1 


1.43 


0.0395 


Gene ontology terms negatively correlated with Luxol fast b\ue-negative white matter area (myelin-loss) 


Generation of precursor metabolites and energy 


5 


DDO, PPOX, GSK3B, ALDH5A1, COXI7 


-1.30 


0.0396 



CDKN2C = cyclin-dependent kinase inhibitor 2C; AFAP1L2 = actin filament associated protein 1 -like 2; ALDH5A1 = aldehyde dehydrogenase 5 family, member A1; 
BLNK = B-cell linker; BTG3 = BTG family, member 3; CALCOC02 = calcium binding and coiled-coil domain 2; CCL2 = chemokine (C-C motif} ligand 2; CCND2 = cyclin D2; 
CD81 =CD81 molecule; CDK2AP1 = CDK2-associated protein 1; CDKN1 A = cyclin-dependent kinase inhibitor 1A; COX1 7 = cytochrome c oxidase, subunit XVII assembly 
protein homolog; DDO = D-aspartate oxidase; DNAJC1 - DnaJ (Hsp40) homolog, subfamily C, member 1; GSK3B = glycogen synthase kinase 3 beta; HBXIP = hepatitis B 
virus x interacting protein; IRAK3 = interleukin-1 receptor-associated kinase 3; KIF23 = kinesin family member 23; LY86 = lymphocyte antigen 86; NEK6 = NIMA (never in 
mitosis gene a)-related kinase 6; PPOX = protoporphyrinogen oxidase; PSMB10 = proteasome (prosome, macropain) subunit, beta type, 10; SMC4 = structural 
maintenance of chromosomes 4; TIMP1 = tissue inhibitor of metalloproteinase 1; TLR1 = toll-like receptor 1; TLR3 = toll-like receptor 3; TOP2A = topoisomerase (DNA) II 
alpha 170 kDa; TPX2 = TPX2, microtubule-associated, homolog. 
doi:1 0.1 371 /journal.pone.009591 7.t004 
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Figure 5. Detection of CDV nucleoprotein by microarray 
analysis, RT-qPCR and immunohistology. (A) Microarray analysis 
showed a low to moderate normalized CDV nucleoprotein RNA 
fluorescence signal in all CDV-infected dogs (group 2-4), which 
significantly exceeded the fluorescence signal measured in control 
dogs (group 1 ). BioB = mean expression of 1 .5 pM Escherichia coli biotin 
synthase (bioB) spiked-in hybridization control RNA; BioC = mean 
expression of 5 pM Escherichia coli biotin biosynthesis protein (bioC) 
spiked-in hybridization control RNA; BioD = mean expression of 25 pM 
Escherichia coli dethiobiotin synthetase (bioD) spiked-in hybridization 
control RNA. (B) RT-qPCR demonstrated robust normalized CDV 
nucleoprotein RNA copy numbers in the cerebella of all CDV-infected 
dogs (group 2-4), whereas the control dogs (group 1) were negative. 
(C) Immunohistology revealed a significantly increased density of CDV 
nucleoprotein antigen-positive cells within the lesions of all CDV- 
infected dogs (group 2-4), whereas the cerebella of the control dogs 
(group 1) lacked a CDV-specific signal. (A-C) Each point represents the 
cerebellar specimen of an individual dog. The median of each group is 
represented by a horizontal line. Significant differences between the 
groups as revealed by the Kruskal-Wallis test with independent pairwise 
post-hoc Mann-Whitney U-tests are marked as follows: *p£0.05. 
doi:1 0.1 371 /journal.pone.009591 7.g005 

by selecting the probe set with the highest significant absolute fold 
change in the respective post hoc test. 

The lists of DEPs or DEGs were checked for significantly 
overrepresented functional terms of the biological process category 



of the gene ontology database employing a modified Fisher exact 
test (EASE score) in DAVID 6.7 [40,41]. Due to the generally low 
frequency of deposited functional gene ontology associations for 
the canine genes, an alternative approach employing orthologous 
mouse genes was used [42] . The resulting lists of gene ontology 
terms were agglomerated into a manageable number of £10 
enriched biological modules employing the DAVID functional 
annotation clustering algorithm [43]. The enrichment score used 
to rank these modules according to their biological relevance is 
calculated as the negative logio of the geometric mean of the 
EASE scores of all incorporated gene ontology terms. 

Gene Set Enrichment Analysis 

For the identification of biological processes intimately associ- 
ated with myelin loss, Gene Set Enrichment Analysis (GSEA, 
Version 2.0.10) was performed employing Pearson's correlation 
coefficient as metric to rank the genes according to their 
correlation to demyelination as observed in Luxol fast blue-cresyl 
violet stained sections and checked for enriched gene ontology 
biological process terms from the Molecular Signatures Database 
(MSigDB) Version 3.1 [44]. The percentage of Luxol fast blue- 
negative white matter in relation to the total white matter on each 
slide was calculated as input. The microarray data set was pre- 
filtered for informative probe sets with a q£0.05 in the multigroup 
LIMMA [45]. Orthologous and unique human gene symbols 
(HUGO), required as feature identifiers by Gene Set Enrichment 
Analysis were retrieved employing MADGene [39] . Gene sets in a 
size range from 4 to 10 genes and positively or negatively 
correlated to demyelination with a p£0.05 were accepted as 
significant and ranked according to the normalized enrichment 
score [44]. 

RT-qPCR 

Reverse transcription quantitative polymerase chain reaction 
(RT-qPCR) was performed from the same batch of RNA that was 
used for the microarray analysis. Total RNA was reverse 
transcribed using the Omniscript kit (Qiagen, Hilden, Germany) 
and random primers (Promega, Mannheim, Germany). Quanti- 
tative PCR and data analysis were performed with the Brilliant 
SYBR-green QPCR Core Reagent Kit (Agilent, Santa Clara, CA), 
the Mx3005P QPCR System (Agilent, Santa Clara, CA), and 
calibration by an external standard curve. Used primer pairs were 
specific for CDV nucleoprotein (sense: 5' GCTCTTGGGTTGCAT- 
GAGTT, antisense: 5 ' GCTGTTTC ACCC ATCTGTTG, ampli- 
con size: 83 base pairs), and glyceraldehyde-3 -phosphate dehydrogenase 
(GAPDH) as house-keeping gene (sense: 5' GTCATCAACGG- 
GAAGTCCATCTC, antisense: 5' AACATACTCAGCACCAG- 
CATCAC, amplicon size: 84 base pairs) [46]. The specificity of 
positive reactions was confirmed by melting point analysis. CDV 
nucleoprotein RNA copy numbers were normalized with a normal- 
ization factor based on GAPDH expression. 

Statistical Analysis 

Analysis of data not otherwise specified was performed using 
IBM SPSS Statistics, Version 20 (IBM Corporation, Armonk, NY, 
USA). Used statistical procedures included the non-parametric 
Kruskal-Wallis test with post-hoc independent pairwise Mann- 
Whitney U-tests for the analysis of the Luxol fast blue-negative 
white matter area, CDV nucleoprotein immunoreactive cells, 
microarray and RT-qPCR data of CDV nucleoprotein RNA, and 
microarray data of manually selected glial cell marker genes. 
Generally, a p-value £0.05 was accepted as significant. Graphs 
were generated employing GrapPad Prism, Version 6.01 (Graph- 
Pad Software, La Jolla, CA, USA). 
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Results 

Histological and Immunohistological changes in CDV 
Leukoencephalitis 

Dogs were categorized into four groups based on the results of 
CDV nucleoprotein immunohistochemistry and histological 
changes as observed in following sections taken from the frozen 
tissue specimens used for RNA extraction. Accordingly, none of 
the control dogs exhibited CDV nucleoprotein immunoreactivity, 
whereas focal to multifocal CDV nucleoprotein-positive lesions 
were detectable in the white matter of the cerebellum and brain 
stem of all CDV-infected dogs (data not shown). The control dogs 
revealed no Luxol fast blue-negative white matter and no gliosis 
and/or inflammatory infiltrates within the cerebellum and brain 
stem (Figure 1A). The focal to multifocal lesions in the dogs with 
acute CDV leukoencephalitis were characterized by slightly 
increased numbers of astrocytes and microglia (Figure IB). 
Additionally, the lesions in the dogs with subacute CDV 
leukoencephalitis exhibited loss of Luxol fast blue-positive white 
matter, indicating demyelination (Figure 1C). The lesions in the 
dogs with chronic CDV leukoencephalitis displayed loss of Luxol 
fast blue-positive white matter and perivascular, lymphohistio- 
plasmacytic infiltrates (Figure ID). 

Top-down Analysis of Major transcriptional changes in 
CDV Leukoencephalitis 

This study was performed in order to gain an assumption-free 
overview of the major transcriptional changes in CDV leukoen- 
cephalitis. An initial principal components analysis of the data sets 
revealed clearly separated clusters for controls and CDV-infected 
dogs, whereas the histologically distinct subgroups of CDV 
leukoencephalitis displayed some overlap, suggesting that CDV- 
infection is the dominant factor influencing the variance (Figure 
SI). One-factorial, multigroup analysis of differential expression 
employing LIMMA revealed 1374 probe sets with q^O.05. 
Sequential filtering employing post-hoc pairwise tests with q£0.05 
and a fold change filter of ^2.0 or < — 2.0 revealed a total of 780 
probe sets corresponding to 442 genes considered as differentially 
expressed in at least one of the six post-hoc comparisons (for a list of 
all 780 differentially expressed probe sets [DEPs], see Table SI). 
Hierarchical clustering of these DEPs revealed four clusters with a 
distinct expression profile (Figure 2). The first cluster included 1 1 1 
DEPs (57 genes; Figure 2, cluster 1, yellow bar) and displayed a 
high up-regulation in all CDV-infected dogs. The second cluster 
consisted of 27 DEPs (16 genes; Figure 2, cluster 2, brown bar) and 
showed a high up-regulation in chronic CDV leukoencephalitis 
only. The third cluster included 495 DEPs (301 genes; Figure 2, 
cluster 3, orange bar) and displayed a moderate up-regulation in 
all CDV-infected dogs. Lastly, the fourth cluster comprised 147 
DEPs (82 genes; Figure 2, cluster 4, green bar) and exhibited a 
moderate down-regulation in all CDV-infected dogs. Only few 
and less informative significantly enriched gene ontology terms 
were retrieved for the two clusters of probe sets highly and 
moderately up-regulated in all CDV-infected dogs employing 
canine probe set IDs as input in the gene ontology analysis (data 
not shown). In contrast, our alternative approach of using the 
orthologous mouse gene symbols revealed many significandy 
enriched gene ontology terms for the three clusters of probe sets 
highly and moderately up-regulated as well as moderately down- 
regulated in all CDV-infected dogs. Functional annotation 
clustering was performed to summarize these gene ontology terms 
into biological modules (Figure 2). The cluster of probe sets 
selectively and highly up-regulated in chronic CDV leukoenceph- 
alitis exhibited no significandy enriched gene ontology terms in the 



automated analysis due to insufficient annotation quality. Manual 
inspection revealed that the majority of the probe sets within this 
cluster (21 out of 27) were related to multiple immunoglobulin-like 
expressed sequence tags. Furthermore this cluster contained the delta 
polypeptide of the cluster of differentiation 3 [CDS) antigen (CD3D). 

The results of all six pairwise post-hoc comparisons are 
summarized in Table 2. Accordingly, large numbers of differen- 
tially expressed genes (DEGs) were detected in the three post-hoc 
tests comparing the three subtypes of CDV leukoencephalitis with 
the control dogs. In contrast, only low to moderate numbers of 
DEGs were detected between the three subtypes of CDV 
leukoencephalitis (for a list of all 442 DEGs, see Table S2). More 
than half of the DEGs in the post-hoc tests comparing the CDV- 
infected- with the control dogs were equally up-regulated in all 
three subgroups of CDV leukoencephalitis and were mosdy 
associated with an activation of the adoptive humoral immune 
response, followed by the cell-mediated immune response 
(Figure 3A). In contrast, the largest subset of down-regulated 
DEGs was observed in subacute CDV leukoencephalitis as 
compared to controls. These down-regulated DEGs were func- 
tionally related to the gene ontology terms "intermediate filament 
bundle assembly" and "regulation of neuron differentiation" 
(Figure 3B). 

Bottom-up Analysis of Transcriptional changes of Glial 
Marker Genes 

The results of an independent analysis of the microarray data of 
manually selected glial cell marker genes employing classical 
statistical methods are shown in Table 3. Accordingly, there was a 
significant down-regulation of myelin basic protein (MBP), myelin- 
associated oligodendrocyte basic protein (MOBP) and proteolipid protein 1 
(PLP1), and a similar trend for 2' ,3' -cyclic nucleotide 3' phosphodi- 
esterase (CMP), myelin associated glycoprotein (MAG) and myelin 
oligodendrocyte glycoprotein (MOG) in subacute CDV leukoencephalitis 
as compared to controls. Additionally, MBP and MOBP were 
down-regulated in subacute CDV leukoencephalitis as compared 
to acute CDV leukoencephalitis. Interestingly, MBP, MOBP and 
PLP1 gene expression revealed a significant reversal to normal 
control levels in chronic CDV leukoencephalitis as compared to 
subacute CDV leukoencephalitis. The oligodendroglial progenitor 
cell (OPC) marker genes chondroitin sulfate proteoglycan 4 (CSPG4, 
synonym: MG2) and platelet-derived growth factor receptor, alpha 
polypeptide (PDGFRA), as well as the Schwann cell marker genes 
myelin protein zero (MP£ synonym: P0) and nerve growth factor receptor 
(jXGFR, synonym: p75) displayed no differential expression. Since 
the expression values of CSPG4, MP/^ and MGFR were clearly 
lower than the mean expression of the smallest 1.5 pM spiked- in 
BioB hybridization control RNA and displayed a constant value in 
all animals, it is assumed that the expression of these genes was 
below the detection sensitivity of the assay. There was a mild but 
significant down-regulation of the Schwann cell myelin-specific 
peripheral myelin protein 22 (PMP22) in acute and subacute CDV 
leukoencephalitis versus controls. The astrocyte marker gene glial 
fibrillary acidic protein (GFAP) displayed an up-regulation in all CDV- 
infected- versus control dogs. 

Exploratory Analysis of Pathways Correlated to 
Demyelination 

In order to identify biological processes intimately related to 
demyelination, Gene Set Enrichment Analysis was performed 
employing Pearson's correlation coefficient as metric to rank the 
genes according to their correlation to myelin-loss as observed in 
Luxol fast blue-cresyl violet stained slides. As already reported 
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above, histology revealed a consistent multifocal loss of Luxol fast 
blue-positive white matter within subacute and chronic CDV 
leukoencephalitis. However, due to the variable proportion of the 
area affected by the multifocal demyelinated lesions as compared 
to the normal appearing white matter area present in each 
individual specimen, the percentage of Luxol fast blue-negative 
white matter area relative to the total white matter area present on 
sections direcdy following the tissue used for RNA extraction was 
calculated as input for Gene Set Enrichment Analysis. Accord- 
ingly, subacute and chronic CDV leukoencephalitis exhibited a 
significandy increased percentage of demyelinated white matter 
area as compared to controls and acute CDV leukoencephalitis 
(Figure 4). Gene Set Enrichment Analysis revealed six gene 
ontology terms positively correlated to demyelination and one 
gene ontology term negatively correlated to demyelination 
(Table 4). A Leading Edge Analysis revealed that the three gene 
ontology terms "viral reproduction", "viral infectious cycle", and 
"viral reproductive process", which were all positively correlated 
with demyelination, can be summarized into a cluster based on the 
shared genes chemokine (C-C motif) ligand 2 (CCL2), CD81 molecule 
(CD81), and topoisomera.se (DMA) II alpha 170kDa (T0P2A) (Figure 
S2). The gene ontology term "generation of precursor metabolites 
and energy", which was negatively correlated to demyelination 
included the down-regulated genes D-aspartate oxidase (DDO), 
protoporphyrinogen oxidase (PPOX), glycogen synthase kinase 3 beta 
(GSK3B), aldehyde dehydrogenase 5 family, member Al (ALDH5A1), 
and cytochrome c oxidase copper chaperone (COX 17). 

Detection of CDV Nucleoprotein Employing Microarray 
Analysis, RT-qPCR and Immunohistochemistry 

The canine genome 2.0 array contains seven probe sets targeting 
different RNAs of the CDV (Table SI). In order to verify the 
microarray results, the data obtained by the probe set targeting the 
CDV nucleoprotein were compared to findings obtained by RT-qPCR 
and immunohistology. Accordingly, the probe set targeting the CDV 
nucleoprotein RNA revealed low to moderate normalized fluorescence 
values in all CDV-infected dogs, which were significantly exceeding 
those of the controls (Figure 5 A). RT-qPCR revealed high normalized 
CDV nucleoprotein RNA copy numbers in all CDV-infected dogs as 
compared to no expression in controls (Figure 5B). In contrast, 
expression of the house-keeping gene GAPDH was detected in all dogs 
and confirmed successful RNA isolation and cDNA synthesis (data 
not shown). As already mentioned above, the quantitative analysis of 
the density of intralesional CD V nucleoprotein-immunoreactive cells 
displayed increased numbers of CDV nucleoprotein-positive cells 
within the lesions of all CDV-infected dogs, whereas the controls were 
consistently negative (Figure 5G). Correlation analysis revealed a 
significant, high, positive correlation of the logio-transformed CDV 
nucleoprotein RNA expression as measured by microarray analysis and 
RT-qPCR (Pearson's correlation coefficient, r = 0.772, p<0.001). 
The CDV nucleoprotein antigen expression as measured by 
immunohistochemistry revealed a significant, moderate to high, 
positive correlation to CDV nucleoprotein RNA expression as 
determined by microarray analysis (Pearson's correlation coefficient, 
r = 0.510, p<0.009) and RT-qPCR (Pearson's correlation coeffi- 
cient, r = 0.698, p<0.001), respectively. 

Discussion 

Top-down Analysis of Major Transcriptional changes in 
CDV Leukoencephalitis 

The transcriptional changes within all CDV leukoencephalitis 
subtypes were dominated by up-regulated genes related to the 
gene ontology terms "immunoglobulin-mediated immune re- 



sponse" and "complement activation, classical pathway". Fur- 
thermore, hierarchical cluster analysis revealed a cluster of DEPs 
highly up-regulated in chronic CDV leukoencephalitis only. The 
majority of these DEPs were specific for multiple immunoglobulin-like 
expressed sequence tags. These findings support previous descriptions 
of an intrathecal immunoglobulin synthesis and deposition in the 
inflammatory phase of CDV leukoencephalitis [14,29]. The 
demonstration of significant amounts of anti-myelin autoantibod- 
ies by Vandevelde et al. (1986) substantiates our hypothesis that 
locally produced autoantibodies and complement may be involved 
in the pathogenesis of demyelination in the chronic inflammatory 
subtype of CDV leukoencephalitis [29], similar to Theiler's 
murine encephalomyelitis [47]. 

Another remarkable group of genes up-regulated in all CDV 
leukoencephalitis subtypes was associated with the gene ontology 
term "positive regulation of T-cell-mediated cytotoxicity". This 
finding is in agreement with the diffuse infiltration of the CNS by 
CD8-positive lymphocytes as observed by immunohistochemistry 
in all subtypes of CDV leukoencephalitis [48,49]. Whether the 
cytotoxic T-cell reaction in CDV leukoencephalitis is solely 
directed against CDV antigens or also spreads to autoantigens in 
later phases is currendy unknown. 

The down-regulation of genes related to the gene ontology 
terms "intermediate filament bundle assembly" and "regulation of 
neuron differentiation" in subacute CDV leukoencephalitis points 
towards a possible involvement of neurons and/or axons in the 
pathogenesis of CDV leukoencephalitis. Notably, infection of 
neurons is thought to be a transient feature of early CDV 
leukoencephalitis [6,18]. Furthermore, the observed down-regu- 
lation of neurofilaments is in agreement with previous studies 
showing progressive axonal pathology and loss beginning in the 
acute to subacute subtype of CDV leukoencephalitis [32]. 
Whether, these changes are a feature of primary axonal 
degeneration (inside-out model) or a response to primary 
demyelination leading to secondary axonal degeneration (out- 
side-in model), needs to be elucidated in further studies [32,50]. 

Bottom-up Analysis of Transcriptional changes of Glial 
Cell Marker Genes 

The observed significant down-regulation of myelin proteins like 
MBP, MOBP and PLP1, and a trend towards lower expression of 
CMP, MAG, and MOG in subacute CDV leukoencephalitis only is 
substantiated by immunohistological studies describing a loss of 
MBP and MAG in experimentally-induced CDV leukoencepha- 
litis starting 16—21 days post infection and prior to the onset of 
perivascular inflammation [51]. Comparably, loss of MBP 
immunoreactivity is described in subacute and chronic lesions 
but not acute plaques in naturally-occuring CDV leukoencepha- 
litis [52]. In vitro studies using CDV-infected mixed brain cell 
cultures indicate that demyelination is either direcdy mediated by 
a restricted infection of oligodendrocytes [23], or indirectiy 
triggered by CDV-induced changes in other glial cells [53]. No 
matter whether the initiating event is direcdy or indirectly CDV- 
induced, in vivo as well as in vitro studies suggest that a metabolic 
impairment of oligodendrocytes leading to down-regulation of 
myelin gene synthesis is the main mechanism of demyelination in 
the initial phase of CDV leukoencephalitis, preceding oligoden- 
drocyte loss by a slow and non-apoptotic "dying back" -mechanism 
[54,55]. The apparent discrepancy between a rebound of myelin 
gene expression to control levels in chronic CDV leukoencephalitis 
despite ongoing demyelination suggests that myelin-loss is 
mediated by other mechanisms in the inflammatory phase of the 
disease, supporting the idea of a biphasic mode of demyelination 
in CDV leukoencephalitis [56]. 
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Exploratory Analysis of Pathways Correlated to 
Demyelination 

Gene Set Enrichment Analysis was applied to identify pathways 
intimately associated with demyelination. Accordingly, the three 
gene ontology terms with the highest positive correlation to 
demyelination were related to gene expression changes known to 
be induced in the host organism by virus replication. This finding 
is supported by the present and previously performed spatiotem- 
poral studies of CDV RNA and protein expression in CDV 
leukoencephalitis, showing that the presence of CDV is a 
prerequisite for the sequential development of demyelination 
[16]. The up-regulated genes shared by all three gene ontology 
terms were CCL2, CD81, and T0P2A, and represent highly 
interesting candidate genes concerning the initial steps of CDV- 
mediated demyelination. CCL2 (synonym: MCP-1) is a pro- 
inflammatory chemokine known to be produced by macrophages, 
T-cells, fibroblasts, keratinocytes, endothelial cells and astrocytes 
that attracts and activates monocytes [57,58]. A toll-like receptor 3 
(TLR3)-dependent up-regulation of CCL2 has been demonstrated 
for Theiler's murine encephalomyelitis virus-infected astrocytes 
in vitro [58]. Furthermore, the observation of increased disease 
severity in Theiler's murine encephalomyelitis virus-infected 
transgenic mice overexpressing CCL2 under the control of a 
GFAP promoter highlights a possible role of virus-induced, 
astrocyte-derived CCL2 in the pathogenesis of demyelination 
[59]. Notably, TLR3, known to recognize dsRNA associated with 
viral infections [5 7] , was also one of the up-regulated genes in the 
present study and is included in the gene ontology term 
"regulation of protein metabolic process". 

The gene ontology term "generation of precursor metabolites 
and energy", which was the only gene set negatively correlated to 
demyelination is defined in the gene ontology database to include 
"chemical reactions and pathways resulting in the formation of 
precursor metabolites, substances from which energy is derived, 
and any process involved in the liberation of energy from these 
substances" [41]. In detail, it comprised the genes D-aspartate 
oxidase (DDO), a peroxisomal flavoprotein that catalyzes the 
oxidative deamination of D-aspartate and N-methyl D-aspartate, 
protoporphyrinogen oxidase (PPOX) involved in heme biosynthesis, the 
serine-threonine kinase glycogen synthase kinase 3 beta (GSK3B) 
participating in energy metabolism and neuronal cell develop- 
ment, the mitochondrial succinic semialdehyde dehydrogenase 
aldehyde dehydrogenase 5 family, member Al (ALDH5A1) that plays a 
role in GABA-metabolism, and cytochrome c oxidase copper chaperone 
(COX17), a non-structural subunit of the mitochondrial respiratory 
chain cytochrome c oxidase [60] . Whether some or all of these 
genes which were collectively down-regulated in CDV leukoen- 
cephalitis are involved in the suggested initial metabolic dysfunc- 
tion of oligodendrocytes needs to be addressed in future studies. 

In summary, the transcriptional changes in CDV leukoenceph- 
alitis support the concept of a biphasic mode of demyelination with 
a presumably initial non-apoptotic oligodendrocyte dystrophy, 
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